Electron-energy bunching in laser-driven soft recoUisions 
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We introduce soft recoUisions in laser-matter interaction. They are characterized by the electron 
missing the ion upon recollision in contrast to the well-known head-on collisions responsible for 
high-harmonic generation or above-threshold ionization. We demonstrate analytically that soft 
recoUisions can cause a bunching of photo-electron energies through which a series of low-energy 
peaks emerges in the electron yield along the laser polarization axis. This peak sequence is universal, 
it does not depend on the binding potential, and is found below an excess energy of one fifth of the 
ponderomotive energy. 
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Recollision of an electron with its parent ion under a lin- 
early polarized strong laser field has been shown to be 
the basis of a plethora of phenomena in atoms [1, 2], 
molecules [3], clusters [4] and solids [5]. In principle the 
recollision process is very simple and a single degree of 
freedom along the laser polarization axis is sufficient to 
describe it (often referred to as the three-step model [1]): 
Firstly, the bound electron is released from an atom due 
to the strong electric field of a laser. Secondly, it is ac- 
celerated and driven back to the ion. In the third step it 
either recombines in the atomic potential or is scattered 
from it. In the former case, high-order harmonics are gen- 
erated (HHG) due to recombination of the electron [6]. 
In the latter case, the elastic head-on collision induces 
the high-energy phenomenon of above-threshold ioniza- 
tion (ATI) with fast electrons emitted [7, 8]. The enor- 
mous impact of HHG up to recent proposals for imaging 
of molecular orbitals [3] and the generation of attosecond 
pulses [9] is not the least due to the simple yet accurate 
description with the three-step model. 

Recently, a surprising strong peak — the "low-energy 
structure" (LES) — was observed at few eV in the photo- 
electron spectrum of atoms in strong infra-red (a few /im 
wavelength) laser pulses [10, 11] and confirmed numer- 
ically with classical calculations [12, 13]. Although the 
LES peak contains about half of the photo electrons it 
was not seen in any of the numerous experiments done 
with 800 nm laser pulses. 

Here, we will give an analytical explanation of the LES 
by introducing a low-energy soft-recollision mechanism. 
It gives rise to a universal series of low-energy peaks in 
the momentum spectrum of the photo electron with well 
defined relative positions of 3/5, 5/7, 7/9 .. . on an abso- 
lute energy scale of about one fifth of the ponderomotive 
energy F^/(4u;^), where F is the amplitude and uj the 
frequency of the laser field. These peaks do not require 
a special binding potential, e. g., long range, nor do they 
need more than one degree of freedom to appear, and 
they can be derived classically since they rely essentially 
on the well known strong- field trajectories as will become 
clear later. 

We will begin by working out the classical structures 



which are responsible for the LES [10], i.e., we consider 
a Hamiltonian H = Hq + V with (throughout the paper 
we use atomic units unless stated otherwise) 

iJo = p2/2 + z F cos{(jjt) (la) 
V = -ll{p'WY^\ (lb) 

describing an electron with position r=(p, z) and momen- 
tum 'P={pp^Pz) using cylindrical coordinates. The elec- 
tron is exposed to the potential V and driven by a laser 
field linearly polarized along z. The probability to mea- 
sure a photo electron with momentum P^, ejected along 
the laser polarization axis z, is given as a two-dimensional 
integral over initial phase-space variables (denoted with 
a prime), 

p(p,) = 11 d^'dpx<^',p;)^(p,-p, p;)) , (2) 

where (j)' = ujt' is the phase of the laser at the time 
when the electron tunnels and is the initial momentum 
perpendicular to the tunneling direction z. The weight 
w{(j)' ^Pp) accounts for the tunnel probability and Jacobian 
factors. The relevant dynamical object in a classical dy- 
namical theory is the deflection function Pz{<i)\Pp) which 
relates final variables to the initial conditions of the tra- 
jectory [14]. Figure 1 shows Pz as a function of the initial 
phase (/)' (or rather the corresponding vector potential 
= Asmcl)' with A = F/u) and the initial transverse 
momentum p^. One can see that Pz develops "finger- like" 
structures with increasing time. They emerge first in the 
second laser period and with each period an additional 
finger appears. These fingers are due to head-on col- 
lisions, responsible for the well-known high-energy phe- 
nomena such as HHG and ATI. Also, these regions are 
characterized by chaotic dynamics [15], very sensitive to 
initial conditions. Responsible for the distinct peaks in 
V{Pz) in Fig. Id at low energies, however, are not the fin- 
gers, but the crossings of contour lines above the fingers 
in Fig. lb, c. They represent saddle points in Pz{4^' ^p'p) 
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FIG. 1: (Color online) The deflection function pz for the momentum along the laser axis as a function of the initial vector 
potential A' —Asincj)' and the initial transversal momentum p'p shown after a) one, b) two, and c) three laser cycles. Spectra 
according to Eq. (2) with w = 1 integrated over the initial parameter range of panels a-c are shown in the right panel d). The 
arrows point to the predicted momenta according to Eq. (7). The laser has a wavelength of A=2 /xm (c(;=0.0228 a.u.) at an 
intensity of /^lO^^W/cm^ (F=0.0534 a.u.). 



with 



dp,/d^' = 0, dp,/dp'^ = 0, 
{d'p,/d^'') {d'p,/dp'') < 0. 



(3) 



Such two-dimensional saddles are known to produce 
peaks since they represent integrable singularities in the 
spectrum [16]. With each additional laser period (for the 
first three we show contour plots and for the first four 
we show the spectra in Fig. 1) a new peak appears. This 
establishes, that the LES actually consists of a series of 
peaks converging towards threshold = 0. 

Figure 2a illustrates that these peaks are due to an 
energy-bunching mechanism of neighboring trajectories: 
The three trajectories shown, start at similar but differ- 
ent phases (j)^ with corresponding drift momenta Asin^^ 
However, they carry the same momentum after the "soft 
recollision" with the nucleus in the second laser period at 
(j)"^ ~ Stt. The trajectory shown in Fig. 2b (corresponding 
to the central trajectory of Fig. 2a) reveals, that it is a 
soft recollision with the nucleus which leads to the energy 
bunching: The electron "misses" the nucleus {p'^ > 0) 
and recollides by virtue of the laser force which turns the 
electron around at z'^ with \z'^\ <^ z = F/u'^ . Hence, this 
new type of a recollision is quite different from the elastic 
reflection off the potential with finite momentum as 
in the head-on collisions in ATI or HHG. How do these 
soft recollisions provide a series of peaks? 

With the characteristics of the soft recollision (all re- 
lated quantities are denoted with a star) as observed 

= z((/)*) - and pI = ((/)*) = , (4) 

this is easy to see using strong-field trajectories 

z{(l)) = z([^'+^^]0 + cos0-l) (5a) 



p^((^) = - Asinc^. 



(5b) 



where we have linearized the solutions of Hamilton's 
equations for Hq in (j)' since tunneling occurs near the 



maximum of the field Fcoscj)' ~ P, i.e., (j)' <C 1. More- 
over, (j)^ = Ap/A accounts for an overall Ap offset of the 
drift momentum due to the potential. As can be seen in 
Fig. 2a this offset does not dependent on ^^ 

From Pz{(l)^) = we get with Eq. (5b) immediately 
= mTT + (-1)'^[0'+^^]. A little thought reveals that 
only odd integers m = 2n + 1 yield non-trivial solutions. 
Requiring that = with the recollision location 



z((/)^) = z[(2n+l)^[(/)'+(/)^] - 2] = 



(6) 



gives the initial phases and in turn the drift momenta 
Pn = 74[(/)^+(/)^], cf. Eq. (5b), for the soft recollisions 



Pn 



F/uj 



(n+l/2)7r ' 



(7) 



From this equation we expect a series of photo-electron 
momentum peaks. This is indeed confirmed by the spec- 
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FIG. 2: (Color online) Time-dependent drift momentum 
Pz{(l>) + Asin((/)) for three trajectories (left panel) with differ- 
ent initial drift momenta showing the effect of the Coulomb 
potential after release [a] and the bunching during the soft rec- 
ollision [c]. Sketch of the rescattering trajectory (right) in a 
Coulomb potential (orange-shaded area). Full trajectory (up- 
per panel) and details (lower) of the three interactions events: 
[a] emission at ~ 0, [b] effectless passing at ~ 37r/2, and 
[cl soft recollison at 6 ^ Stv. 



tra shown in Fig. Id, where peaks appear cycle by cycle. 
How do these peaks emerge? 

From our analysis so far, one-dimensional (ID) dy- 
namics with some potential (short or long range) should 
be sufficient to explain the underlying mechanism. To 
this end we consider the ID Hamiltonian H = j9^/2 + 
xF cos{ujt) + Vs{x)^ with position x, momentum p and 
the range s of the potential. Starting always at the ori- 
gin X = 0, but with different phases the electron 
is propagated until \x\ s and the drift momentum 
p{(j))-\-Asm{(j)) is constant. The deflection function p{(l)') 
along with the corresponding photo-electron spectrum 
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(8a) 
(8b) 
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is shown in Fig. 3 for the short-range Gaussian potential 

V;(x) = -exp(-(x/s)2)/s. (9) 

By integrating the force due to the potential Vq the de- 
flection function can be written in the form 



1 /"^ 
W) = — / 



Ax 



(10a) 
(10b) 



where A(j)' (shown as a dashed line in Fig. 3a) represents 
the contribution from the laser field without the poten- 
tial. The second term 5p represents the impact from the 
external potential which leads to modulations in p{(t)'). 
Whether the modulations are really visible as pronounced 
peaks in the spectrum depends on the strength of the im- 
pact 5p{(l)') in Eq. (10b). 

Peaks occur in the first place if dp/d(j)' = 0, cf. (8b). 
Physically, this means that the change in the impact 
strength d5p/d(j)' must exactly compensate the change 
in the drift momemtum, which is simply A. A weak im- 
pact leads only to a marginal decrease of the slope of 
p{(j)') giving rise to a shallow hump in the spectrum. On 
the other hand a strong impact 5p will overcompensate 
the change of the drift momentum leading to a negative 
slope for some (j)' accompanied by two extrema. This is 
indeed the case in Fig. 3 for the higher recollisions n > 1. 
The cross-over between weak and strong impacts is de- 
termined by d^p/d(j)''^ = 0, which we take as a measure 
for a potential to produce pronounced peaks. These two 
conditions allow us to determine the initial phase 0' and 
the strength parameter s of the potential for producing 
pronounced peaks through soft collisions as we will show 
now analytically. To this end we consider the integral, 
cf. Eq. (10b) 
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FIG. 3: (Color online) Deflection function (left) and cor- 
responding photo-electron spectrum with finite resolution 
(right) for the Gaussian potential Eq. (9), laser parameter as 
in Fig. 1. The missing interval in the deflection function rep- 
resents initial conditions which lead to trapped trajectories. 
The dashed line corresponds to the strong field drift momen- 
tum p{(j)') = A sin (/)' . 



for a soft-recollision trajectory 



(lib) 



It is defined by the quadratic dependence of a strong- 
field trajectory around the recollision phase 0*. One can 
extend this behaviour to ^ ^ ±oo, since the force in 
Eq. (11a) vanishes for large Note that the impact Sp 
depends through the recollision point x'^ on the initial 
phase (j)' . 

Fullfilling the condition d'^p/dcj)^'^ = can be cast into 
the form /2(c^) = for the ratio = x'^/s and the 
function 



kj 



+00 d^+^V,{x) 

CO 



X = C*+(/)2/2 



kj = [(2n+l)7r]^ x- 



(12a) 
(12b) 



This follows directly from Eq. (11) by using the chain rule 
and an appropriate rescaling of the integration variable. 
Note, that the integral in Eq. (12a) does only depend on 
the shape of the potential, but not on any specific pa- 
rameters of the problem 5, F, or uj. Hence, the value 
which solves f2{c^) = is a general constant which 
assumes the value = —0.319 for the Gaussian poten- 
tial (9). The first condition dp/d(j)' = reads with the 
definition Eq. (12) simply A + fi{c^) = which can be 
recast into a form that determines 5 as a function of the 
laser parameters F, u and the order of the recollision n. 



[2(2n+l)7r/i(c* 
(Fw2)i/5 



12/5 



(13) 



X = x{(f)) 



This allows us to determine quantitatively the scale s and 
through the relation x^ / s = also the point of the recol- 
lision x'^ at which the deflection function has a zero-slope 
inflection point. Our quasi-analytical determination of 
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TABLE I: Soft-recollision parameters (in a.u.) for the Gaus- 
sian (9) and the soft-core (14) potential as obtained from 
Eqs. (13) and (16), respectively. For comparison full numeri- 
cal results for propagation until t = 2T and t = 3T, respec- 
tively, are shown in italics for the laser parameters of Fig. 1. 



ID Gaussian 


ID soft 


-core 


3D Coulomb 


n = l 


n = 2 


n = 1 


n = 2 




n = l 


s 32.7 


40.2 


24.6 


30.2 




23.6 


31.9 


38.4 


23.9 


29.4 




22.2 


-X* 10.5 


12.8 


6.5 


8.0 




10.9 


10.2 


12.6 


6.4 


8.1 




12.4 



the soft-collision parameters s and x'^ is remarkably ac- 
curate as the comparison with the numerical exact values 
from the soft colliding trajectory propagated under H 
reveals in Table I. There, we also list the corresponding 
values for the ID soft-core Coulomb potential 



Vs{x) = -l/ {x'^+s') 



1/2 



(14) 



for which one obtains = —0.264. 

In fact, the 3D physical case discussed in the beginning 
can be mapped onto the ID soft-core potential since p is 
very slowly varying across the soft collision (see Fig. 2b) 
and can be effectively treated as a parameter, i.e., we 
take at the soft collision = s m. the soft-core potential 
(14). In order to fullfill the saddle-point conditions (3) 
we only have to exchange d?p/d(j)''^ = from our ID 
treatment with dpz/dp = 0. The latter reads 



/ + 00 
-oo 



2— z — 
dz dz-^ 



= 0, 



(15) 



2 = C*+</)2/2 



and can be expressed with integrals from Eq. (12) produc- 
ing an equation only dependent on = z^/p*. A sim- 
ilar procedure as described above for the ID case yields 



= -0.462 and ultimately 
p* = 2.90/ {Fco'Y^' 



(16) 



for the first recoUision (n = 1) in very good agreement 
with the numerical values, see Table 1. 

In summary, we have identified a soft-recollision mech- 
anism which induces energy bunching for low- energy 
photo electrons along the laser polarization. The bunch- 
ing occurs since electrons with initially different drift mo- 
menta can aquire impacts through soft recollisions which 
exactly counterbalance the initial differences leading to 
a series of photo-electron peaks with relative positions 
Pn/Pn-i = (2n— l)/(2n+l) for n>l. This is a universal 
result which does neither depend on the dimensionality 
of the potential (one degree of freedom is enough) nor on 
the character of the potential (short or long range) or the 
laser intensity and frequency. It does, however, require 
a quiver amplitude z much larger than the range of the 
potential s. This is necessary to provide well defined im- 
pacts Sp by the potential when the mainly laser driven 
electron trajectory passes the potential. 

The absolute positions Pn of the peaks are slightly de- 
pendent on the potential and the laser pulse and will also 
be infiuenced by focal- volume averaging and the pulse 
envelope. This applies in particular to the higher-order 
peaks very close to threshold where also additional dy- 
namical effects may mask the soft recollision peaks. How- 
ever, since each peak pn is generated in a successively 
later laser period n+1, one can in principle control the 
number of peaks by varying the total number of cycles 
in the laser pulse [17]. The phenomenon is essentially 
classical because the potential perturbs the strong field 
dynamics only marginally. The latter contains only up to 
quadratic operators (depending on the length or veloc- 
ity gauge) . Owing to the Ehrenfest theorem the quantum 
evolution can therefore be described equivalently by clas- 
sical mechanics. 



[1] P. B. Corkum, Phys. Rev. Lett. 71, 1994 (1993). 

[2] M. Y. Kuchiev, JETP Lett. 45, 404 (1987). 

[3] J. Itatani et al, Nature 432, 867 (2004). 

[4] U. Saalmann and J. M. Rost, Phys. Rev. Lett. 100, 

133006 (2008). 
[5] F. Brunei, Phys. Rev. Lett 59, 52 (1987). 
[6] M. Lewenstein, P. Balcou, M. Y. Ivanov, A. L'Huillier, 

and P. B. Corkum, Phys. Rev. A 49, 2117 (1994). 
[7] G. G. Paulus, W. Becker, W. Nicklich, and H. Walther, 

J. Phys. B 27, L703 (1994). 
[8] D. B. Milosevic, G. G. Paulus, D. Bauer, and W. Becker, 

J. Phys. B 39, R203 (2006). 
[9] F. Krausz and M. Ivanov, Rev. Mod. Phys. 81, 163 

(2009). 



[10] C. I. Blaga et a/., Nat. Phys. 5, 335 (2009). 
[11] W. Quan et al, Phys. Rev. Lett. 103, 093001 (2009). 
[12] C. Liu and K. Z. Hatsagortsyan, Phys. Rev. Lett. 105, 
113003 (2010). 

[13] T.-M. Yan, S. V. Popruzhenko, M. J. J. Vrakking, and 
D. Bauer, Phys. Rev. Lett. 105, 253002 (2010). 

[14] J. M. Rost, Phys. Rep. 297, 271 (1998). 

[15] G. van de Sand and J. M. Rost, Phys. Rev. Lett. 83, 524 
(1999). 

[16] The situation is similar to the logarithmic singularities of 
the density of states in a two-dimensional crystal: L. V. 
Hove, Phys. Rev. 89, 1189 (1953). 

[17] A. Kastner, U. Saalmann, and J. M. Rost, in preparation 



